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I.   Introduction 

This  report  describes  a  two  dimensional  stochastic  process 
and  a  program  for  its  simulation.   The  process  can  be  considered 
to  describe  a  random  track  of  predetermined  length  that  origi- 
nates and  terminates  at  predetermined  end  points.   The  process 
is  a  modification  of  a  process  that  is  without  the  length  and 
termination  constraints.   This  process,  which  is  also  described, 
is  a  generalization  of  the  random  tour  process  described  by 
Washburn  in  Reference  1. 

A  particular  process  is  defined  by  the  track  length,  the 
maximum  expected  track  segment  length  and  the  distance  between 
the  end  points.   For  motion  at  a  constant  speed,  the  track 
length  constraint  is  equivalent  to  a  transit  time  constraint. 

In  Section  II,  the  process  and  the  unmodified  process  are 
described  and  some  of  their  characteristics  are  discussed.   In 
Section  III,  the  simulation  program  is  described.   By  using 
the  program,  graphs  and  data  can  be  produced  that  represent 
realizations  of  the  process.   The  program,  which  is  written  in 
BASIC,  is  for  an  HP-85. 


II.   A  Description  of  the  Process 

The  process,  as  it  is  described  here,  generates  a  random 
track  of  predetermined  length  that  connects  two  predetermined 
end  points.   The  track  originates  at  a  predetermined  starting 
point,  the  first  point,  and  terminates  at  a  predetermined  stopping 
point,  the  final  point.   The  track  is  continuous  between  these 
two  points  and  consists  of  a  sequence  of  straight  line  track 
segments.   Some  representative  tracks  are  shown  in  Figure  1. 

A  sequence  of  rectangular  coordinate  systems  is  used  to 
define  the  process.   The  first  is  a  reference  system  for  the 
first  track  segment,  the  second  is  a  reference  system  for  the 
second  track  segment  and  so  on.   The  x-axis  of  the  first  system 
is  coincident  with  the  line  joining  the  first  and  final  points, 
the  origin  is  at  the  mid-point  of  the  line  and  the  system  is 
oriented  so  that  the  final  point  is  on  its  positive  x-axis. 
For  each  successive  system,  the  x-axis  is  coincident  with  the 
line  joining  the  last  point  of  the  preceding  track  segment  and 
the  final  point,  the  origin  is  at  the  mid-point  of  the  line 
and  the  system  is  oriented  so  that  the  final  point  is  on  its 
positive  x-axis.   Each  track  segment  is  determined  by  two  quan- 
tities:  an  angle  9  and  a  length  r.   The  angle  6  is  measured 
clockwise  from  the  direction  of  the  positive  y-axis  of  the 
reference  coordinate  system  of  the  track  segment.   It  is  the 
relative  course  that  will  generate  the  track  segment.   The 
length  r  determines  the  track  segment  length. 
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'igure  1.   Some  representative  trades.   The  tracks  were  generated 
by  using  the  program  that  is  described  in  Section  III. 


The  process  is  defined  as  follows:   For  each  track  segment 
except  the  last  a  relative  course  9  is  chosen  from  a  distribu- 
tion with  density  function 

f0(8)  =  (1/2tt)  b/(a  -  c  sin  6) 

where  0  <  9  <  2tt,  2c  is  the  distance  between  the  first  and 

final  points,  2a  is  the  track  length  with  a  >  c    0  and 

2     2    k 

b  =  (a   -  c  )  .   Next,  a  length  r  is  chosen  from  a  conditional 

exponential  distribution  with  density  function 


fR>Q(r|9)  =  [l/6(9)]exp[-  r/5(9)] 

where  0  <  r,  6(9)  =  (a  -  c)6/(a  -  c  sin  9)  and  -5  is  the  maximum 

expected  track  segment  length.   Using  the  parameters  a,  ,  b,  and 

2 

c,  defined  below,  if  r  >  2a,  ,  then  r,  =  b,  /(a,  -c,  sin  8,  )  .   In 

this  case,  the  (k+l)st  track  segment  is  the  line  that  joins   the 

last  point  of  the  kth  track  segment  to  the  final  point  and  the 

orocess  terminates.   If  r  <  2a,  ,  then  r,  =  r.   In  this  case,  the 

k        k 

last   point  of  the  kth  track  segment  and  the  final  point  deter- 
mine the  (k+l)st  reference  coordinate  system  and  the  process 
continues.   The  parameters  in  the  above  expressions  are  defined 
as  follows:   2a,  is  the  remaining  track  length  with  2a,  =  2a 

and  2a,  =  2a,  , -r,  ,  for  k  >  1 ;  2c,  is  the  distance  from  the 
k     k-1   k-1  k 

last  point  of  the  (k-l)st  track  segment  to  the  final  point  with 

2    2  3" 

2c,  =  2c  and  b,  =  (a,  -c,  )  2.   The  program  that  is  described  in 

j_  k      k    k  3 

Section  III  can  be  used  to  simulate  this  process. 

Figure  2  illustrates  the  generation  of  a  three  segment 
track.   The  ellipse  on  which  the  second  segment  terminates 


First 
Point 


Final 
Point 


Figure  2.   A  three  segment  track  and  a  focusing  ellipse 
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Figure  3.   The  first  rectangular  coordinate  system,  the 
associated  polar  coordinate  system  and  the 
reference  ellipse. 


and  the  third  segment  begins  has  one  focus  at  the  last  point  of 
the  first  segment  and  the  other  at  the  final  point.   With  2c„ 

the  distance  between  these  two  points,  2a~  =  2a  -  r,  the 

2     2      2 
remaining  track  length  and  b-   =  ao   ~  co  '  ^e  ellipse  is 

2    2     2    2 

defined  through   x  /a~   +  y  /b0   =  1  in  the  second  reference 

2 

system.   It  is  also  defined  by  r  =  b-  /(a~  -  c~  sin  9)  in  the 

polar  coordinate  system  with  its  origin  at  the  last  point  of 
the  first  segment  and  its  polar  axis  parallel  to  the  y-axis  of 
the  second  reference  system. 

As  a  basis  for  discussing  the  relative  course  and  track 
length  distributions  that  are  associated  with  the  process, 
consider  the  polar  coordinate  system  with  its  origin  at  the 
first  point  and  its  polar  axis  parallel  to  the  y-axis  of  the 
first  reference  system.   Next,  consider  a  reference  ellipse 
with  one  focus  at  the  first  point  and  the  other  at  the  final 
point  that  is  defined  by  r  =  b  /(a  -  c  sin  9) .   The  coordinate 
system  and  the  ellipse  are  illustrated  in  Figure  3. 

The  distribution  of  the  relative  course  3  and  the  condi- 
tional exponential  distribution  of  the  length  R  are  both  deter- 
mined by  p(9)  =  (a-c)/(a  -  c  sin  9) ,  the  ratio  of  the  value  of 
the  radial  coordinate  of  the  ellipse  at  9  to  its  maximum  value 
at  tt/2.   Note  that  when  a  track  segment's  relative  course  is 
tt/2  the  segment  is  directed  toward  the  final  point.   The  dis- 
tribution of  the  relative  course  0  is  determined  by  the  require- 
ment that  the  ratio  of  the  probability  that  a  relative  course 
is  between  9  and  9  +  A8  to  A9  is  proportional  to  o(9) .   The 


conditional  exponential  distribution  of  the  length  R  is  deter- 
mined by  the  requirement  that  the  ratio  of  5(6),  the  expected 
track  segirent  length  for  9  ,  to  5 ,  the  expected  track  segment 
length  for  it/2,  equals  p(9). 

If  the  first  point  and  final  point  are  coincident,  then 
c  =  0  and  the  reference  ellipse  becomes  a  circle.   In  this  case, 
the  distribution  of  the  relative  course  0  is  uniform  with  den- 
sity function  fQ(9)  =  1/2tt  where  0  <  9  <  2rr,  the  relative  course 
9  and  the  length  R  are  independent  and  the  distribution  of 
R  is  exponential  with  density  function  fp(r)  =  (1/5) exp  (-r/5) 
where  0  <  r. 

Consider  a  process  that  is  defined  in  the  same  way  except 
that  the  length  of  a  track  segment  is  always  chosen  to  be  equal 
to  the  value  determined  for  R.   For  this  process,  a  is  not  the 
track  length,  there  is  not  a  length  constraint,  and  the  final 
point  is  not  an  end  point,  there  is  not  a  termination  constraint. 
The  final  point  is  a  reference  point  of  the  process,  literally,  a 
focal   point  of  the  process.   If  c  =  0,  the  process  is  identical 
to  the  random  tour  process  described  by  Washburn  in  Reference  1. 
If  c  >  0,  it  is  a  generalization  of  that  process. 

By  inspection  of  the  density  function,  it  can  be  seen  that 
both  the  mean  and  the  mode  of  the  relative  course  ?  are  equal 
to  tt/2.   For  the  constrained  process,  this  implies  that  the 
most  likely  next  to  last  and  last  track  segment  relative  course 
values  are  tt/2  and  -tt/2.   In  Appendix  1,  it  is  shown  that  the 
expected  value  of  the  length  R  is  a6/(a+c).   For  the  uncon- 
strained process,  this  implies  that  the  average  track  segment 
length  is  aS/(a+c). 


A  new  process  could  be  generated  by  redefining  the  para- 
meters a  and  c.   For  example,  one  could  use  the  definition 
a  =  2c  where  c  is  again  the  distance  between  the  first  and  final 
points.   In  this  case  the  degree  of  focusing  of  the  process 
would  be  constant.   However,  relative  to  the  process  described 
here,  it  is  conjectured  that,  for  track  lengths  greater  than  2a, 
the  process  would  concentrate  track  points  about  the  final  point, 
for  track  lengths  less  than  2a,  the  process  would  concentrate 
track  points  about  the  first  point.   In  the  simulation  program, 
for  a  track  of  m  segments  where  m  >  2,  the  maximum  length  of 
the  first  m-2  track  segments  is  computed.   And,  it  is  also  con- 
jectured that,  for  a  track  length  less  than  2a,  the  probability 
that  the  length  of  the  last  track  segment  would  exceed  this 
statistic  would  be  greater  than  that  for  the  process  described 
here. 


III.   The  Simulation  Program 

The  simulation  program  requires  three  inputs:   The  distance 
D  between  the  end  points,  the  track  length  TL  and  the  maximum 
expected  track  segment  length  6.   The  program  is  for  an  HP-85 
and  is  written  in  BASIC  for  that:  computer.   The  program  listing 
is  in  Appendix  2.   After  loading  the  program  and  initiating  RUN, 
a  user  is  prompted  for  these  inputs  and  advised  cf  their  limits. 
An  example  of  the  prompting  with  a  set  of  inputs  is  shown  below. 


DISTANCE  Oi.Di.lQO? 

19 

TRACK  LENGTH  <D<TL158 ) 9 

58 

DELTA  <!_<£>? 


Table  1.   The  input  limits  are  based  on  display  and  memory 
requirements . 

As  a  simulation  progresses,  the  developing  track  is  displayed 
on  the  HP-35  CRT.   After  a  simulation  is  completed,  a  user  is 
given  the  option  of  displaying  track  data  or  printing  the  track 
and  the  track  data.   In  the  CRT  display  of  a  -rack,  the  stalling 
point  is  on  the  left  side  and  the  stopping  point  is  on  the  right 
side.   In  the  printed  output,  a  track  is  rotated  90°  so  that  the 
starting  point  is  at  the  top.   Figure  4  shows  sections  of  the 
printed  output.   In  the  columned  data,  I  is  the  track  segment 
number,  <j>  is  the  track  segment  course  for  a  starting  point  that 
is  due  west  of  the  stopping  point,  R  is  the  track  segment  length, 


9  is  the  track  segment  relative  course  and  X  and  Y  are  the 
coordinates  of  the  last  point  of  the  track  segment  in  the  first 
coordinate  system.   If  the  starting  point  is  considered  to  be 
due  west  of  the  stopping  point,  then  the  positive  x-axis  of  the 
first  coordinate  system  is  directed  east  and  the  positive  y-axis 
is  directed  north.   Recall  the  origin  of  the  first  coordinate  system 
is  at  the  mid-point  of  the  line  joining  the  starting  and  stopping 
points. 

For  the  track  shown  in  Figure  4,  the  termination  criterion 
was  satisfied  at  the  simulation  of  the  49th  track  segment.  In 
the  first  line  below  the  columned  data,  0  is  the  average  track 
segment  relative  course  where  -90 °<  8  270 °,  R  is  the  average 
and  RM  is  the  maximum  track  segment  length  for  the  first  49 
track  segments.  In  the  second  line,  $(51)  is  the  relative  course 
of  the  51st  and  last  track  segment  and  R(51)  is  its  length. 

If  a  track  simulation  reaches  97  segments  without  termina- 
tion and  the  termination  criterion  is  not  satisfied  for  the  98th 
segment,  the  simulation  is  stopped.   In  this  case,  the  values 
on  the  first  line  below  the  columned  data  refer  to  all  of  the 
98  track  segments  and  "TRUNCATED"  is  on  the  second  line.   For 
both  cases,  the  average  relative  course  9,  the  average  track 
length  R  and  the  maximum  track  length  RM  relate  to  the  uncon- 
strained stochastic  process. 

Figures  5,  6  and  7  show  sections  of  the  printed  output  for 
three  additional  simulated  tracks.   The  plots  illustrate  the 
effect  of  the  value  of  5  on  the  plot  scale.   Doubling  the  value 
of  o  doubles  the  distance  between  points  on  the  plot.   This 


10 


increases  the  likelihood  that  the  HP-85  CRT  will  display  a  com- 
plete track.   The  plots  illustrate,  to  some  extent,  the  increased 
focusing  that  occurs  as  the  value  of  TL  approaches  the  value  of 
D. 

A  second  simulation  program  to  generate  track  data  in  a  non- 
graphical  but  more  efficient  manner  was  written.   A  preliminary 
analysis  of  some  data  generated  by  the  simulation  has  been  made. 
In  Table  2  below,  some  estimates  p  with  associated  95%  confidence 
intervals  are  given  for  p,  the  probability  that  the  last  track 
segment  length  will  exceed  the  maximum  track  segment  length  of 
the  preceding  track  segments. 


TL  6  p  95%  confidence 

interval 


10  50  1  .60  .57  <  p  <  .63 

10  50  2  .55  . 52  <  p  <  .58 

10  50  5  .53  . 50  <  p  <  .57 

10  20  1  .44  .40  <  p  <  .47 

10  20  2  .43  . 39  <  p  <  .46 

10  20  5  .45  . 42  p  <  .48 

Table  2.   Estimates  p  for  p,  the  probability  that  the  last 

track  segment  length  will  exceed  the  maximum  track 
segment  length  of  the  preceding  track  segments. 
Each  estimate  is  based  on  a  sample  size  of  900. 
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Figure  4.   A  program  generated  track.   The  track  data  for 
segments  4  through  4  7  have  been  omitted. 
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Figure    5.      A   program   generated    track. 
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Figure    6.      A   program  generated    track 
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Figure    7.      A  program  generated   track 


Appendix  1   Some  Mathematical  Relationships 

In  the  simulation,  values  of  0  are  generated  by  using  the 

inverse  transform  method.   This  method  requires  the  use  of 

the  cumulative  distribution  function  Ffl(8)  and  it  can  be  found 

as    follows:       For      -tt/2    <    9    <    tt/2    , 

9 
FQ(6)    =      /    (1/2tt)  [b/(a-c    sin    9)]    d9 

y  -tt/2 

=  2[b/(a+c) ]{tan"1[ (a  tan  (9/2)-c)/b]  -  tan-1[- (a+c) /b] } 
For  tt/2  <  9  <  3rr/2,  FQ(9)  =  1  -  Fq(tt-8). 

The  expected  value  of  R  for  the  unmodified  process  that  is 

given  in  Section  II  can  be  found  as  follows: 

tt/2 
E(R)     =    2    /      E(R|0=9)     f0(9)     d9 
-tt/2 

tt/2 

/       (5/ir)    b(a-c)  [l/(a-c    sin9)2]    d9 
-tt/2 


=  (25/tt)  [a/ (a+c)  ]{tan  1[b/(a+c)]  +  tan  X[(a+c)/b]} 


=  a  5 /a+c 


If  the  termination  criterion  is  satisfied  on  the  first 

track  segment,  the  track  will  consist  of  only  two  track  segments 

2 
With  r  (9)  =  b  /(a-c  sin  6),  the  probability  P0  of  this  event 
e  £ 

can  be  found  as  follows: 
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TT/2 

P2   =    2      /    P(R    >    r    (8) |0   = 
-tt/2 


8)     fQ(6)    d9 


tt/2 
=   2     /    exp[-re(8)/5(8)]    fQ(8)    df 
-tt/2 

=   exp[-  (a+c) /5] 


Appendix  2.   The  HP-85  Simulation  Program 


1  8 

29 

38 
48 

58 

68 
78 

S  8 

98 
1  88 
118 
128 

138 


148 

158 
168 
178 
188 

1  '38 

2  8  8 
218 
228 


238 


■^ 

48 

£ 

58 

£, 

68 

-■} 

78 

^ 

3  8 

w 

ft  8 

"7 

88 

"7 

18 

"7 

28 

"7 

38 

T 

4y 

- 

58 

7 

68 

^ 

7  8 

T 

58 

"T 

98 

488 

4. 

18 

4 

28 

4 

7  8 

a 

48 

4 

=:£1 

a 

t'K} 

COM  r  '■  1 8  8  )  .■  R  ■■  1 8  8  .'■  ,  T  '•'.  1  8  8  . >  >  '<  <  '■.  1 

8  8  >  •''  '■  1  8  8  '•• 

0 1 3  P  "  D 1 3  T  h  N  C  E  "  ;  "  i  8  "  .;  C  H  R  *  v  1 

33  >  >   "  D  "  ;  CHR $  <  1 33  >  .:  "  1 8  >  "  ; 

INPUT  DS  Cl=C---'2 

IP  0>18  OR  D<8  THEN  28 

3 1 SP  "TRACK  LENGTH  "  ,  "  <  D  "  .:  CH 

RS<68>;  UTL"  ;CHR*U88>i  "58>"  , 

INPUT  L 

IF    L>58    OR    L    =0    THEN    58    ELSE 

ftl=L-'2 
DISP     "DELTA";"     <Li"iCHft*C188> 
.;  CHRKIS);  ">"  ; 
INPUT    Di 

IF    01 <=8    THEN    38 
RANDOMIZE 
v  ,-;  gj ;.,  =  -  Q  i    i»    v  <  g  >  =  y    i*    F  •:.  8  >  =  8    '2 

p  (  8  .» =  8    8    T  <  8  >  =  8 
Bl=8    8    82-0    8    83=8    8    Rl=8    8 
fi=Rl    8    C=C1    8    3=8    8    V=8    8    W— 
8 

GCLEftR  8  E= 1 8£SQR ( 0 1 : 
SCRLE  -E ..£..-<  .  75*E>  •■  .  75*E 
MOVE  CI, 8  8  L ft BEL  "o" 
M  0  V  E  -  C 1  ,  8  8  L  ft  6  E  L  "  o  " 
FOR  1=1  TQ  98 
U-RNO  8  Gl=8 
IF  U  :,=  .  5  THEN  228 
U=U- .5  8  Gl =1 

Nl*SQRt<fll+Cl>^<fll-Cl  J  ^    8  i-iZ 
=  3QR':  i-Ci#Cl-'<Al*fti  )  ) 
T=2*ATN<N2*TAN<PI*U-ATN»:  Nl  >  > 

+  "'  1    •■ ' Hi    ) 

IF    Gl -t     THEN    T=PI-T 

P=  i.  ft  l  -C  i  ••  •■'  i  ft  1  -C  1  *S  I H  (  l)  ) 

02=D1*P 

R=-<  02*LOG<  i-RHCO  ) 

F-T+S 

GOSUB    388 

I r    ft2> C2    THEN    348 

R=<A*A-C#C>'-  f  ft-C*SIH<T>  > 

GOSUB    388 

D  R  ft  W    ',-•', '  I  ."'  •  V '.  1  ) 

ft  =  ft  2    8    C  =  C  2    8    3  =  -  ft  T  N  2  ■-.  V  .■  X  >     5 

H=T  8  GOSUB  778 
T<I>=H 

H=r  8  GOSUB  778 

r  •.  a  -  ■  -  n 

IF  B2=l  THEN  458 

Rl=MftJ<{<R,  Rl  ) 

y  =  V+<T-V>/'I  S  W-W+^R-itf  :'>"I 

NEXT  I 

0I3P  8  DISP  "TRUNCATED"  8  K> 

I  8  B  3  = 1  8  GOSUB  358  8  GOT  0 

438 

0  ft  ft  Ui  C  1  .<  8 

K  =  1  +  1  8  P  '  -;"  >  =  2  3-  C  !*  T ».'  k  ">  =  ft  ft.     * 


!=S+PI--*2  8  GOSUB 


I*   P  ;  k 


479 

489 

4 '?  9 

599 

519 

529 

53W 


GOSUB  350 

DISF  «2  DISP  "DISP  DATR  i  ■•  -  H  ■ 

!1 

INPUT  hi 

IF  R*="Y"  THEN  529 

GOTO  €39 

DISP  USING  539  ;  CHR$<.  15>,CH 

R  %  i  1  6  '? 

T  *rt  Cl  i—  C  T    ••      '-.-'       ii  T  ;i       "7  v      £■       .,  ■ .'       >i  ri  ii 

i  ■•  .  M  »j  t       >  -i  •     .■  i  •.  .•        i        ■  -i  rt  .•  H  .■  ■+ 1\  ..       K 


5^y 
559 

569 

5  79 

5  S  9 

5'?  9 

699 

5  i  9 

6  2  9 


*-, 

: -9 

f^ 

4  9 

6 

59 

K 

69 

6 

79 

s 

S  9 

-, 

39 

"? 

99 

X 

1  9 

■? 

29 

— 

7  i~i 

FOR  1=1  TO  K-l 

DISP  USING  569 

,  T  <  I  .>  ,  A  CI  >  ,  Y  •:  I  > 

IMAGE  ...••-..  2D. -  2X.-3Z..  X..  2D  2D.2X 

,  3Z,;<,  3D  2D,X.3D.  2D 

NEXT  I 

DISP  USING  539  .;  CHR$-.  1  44  .:>  ,  V 

,  CHR*<  219.)  ,  W,  Rl 

I  URGE  ,■■■■'■  R,  "  =  "  ,  3Z  •  2S,  R,  "  ='"  . 

20  2C  .•  2>' ,  "  RM=  "  ..  2D  2D 

Ic  B3-~l  THEN  649 

DISP  USING  629  ?     CHR-f  •:.  1  5  >  ..  K  '  = 

F  <  K  >  •  K ,  R  ;  K  ':■■ 

I  MAGE  ■  ■••- ,  R ,  "  '■  !S  2D ,  "  >  =  "  3Z  >  2X 

..  "  R !-  "  ..  2D  .•  "  >  =  "  >  2D  .20  .■  3--'  ■ 

G 0  T  0  669 

DISP  USING  659  .;  "TRUNCRTED"' 

I*1RGE    ,     ,  9R.<  ?••'  .■ 

I ~    B 1 =  9    THE  N    6  3 9 

CRT     I S     1     ~±    G 0 TO    71 9 


H 


7  49 
7  5  9 
7  69 
77  9 
7  3  9 
799 
3  9  9 
3  i  9 
329 
3  3  9 
349 
359 
ft  f,  9 


DISP  i2  DISP  "PRINT  PLOT  ( 

Tf.  cr.--  n:<  " .; 

I  -••  P  U  T  B  -f 

IP  B*="Y"  THEN  729 

RLPHR  ■'£    END 

PRINT  USING  739    D  ,  L  ..  D  1 

IMAGE  ■  3-  ,  "DISTANCES"  ,  DD  DO  ■ 

,  "TRACK  LENGTHS"  ■  ODD  D  D  .■■••' ••- 
.  ;DELTR=^  .  DD  DO.  3-- 
C R T  IS  2  i?  61  =  1 
GRAPH  i    COP'-' 
G 0 TO  52 9 

?  T  D  '■  H  .'.'■  S  H  =  R  flO ,:  H  ,  369,' 
If     h -;:  9  THEN  H=H-*-369 
RETURN 

><>..  I  ■='A'.   I  -1  >+R*SIN«  F) 
V  >'.   I  ;  =  Y  i.   I  —  I  WPtC 0 S  •  r 
;:<  =  f:  i  —  '■£  (  I  )  !*  V  =  —  Y  •'•'  I  ,'! 
C 2  =  S Q R  '  •- £ X  +  v  '*"  Y  '  / 2  '2  R2  =  R-F  ••"'  I 
RETURN 

H = ',,•  ii  H = R  T  0  '■  H  •'  !2  v = h 
RETURN 
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